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ABSTRACT 

We present a multi-wavelength study of the emission-line nebulae located ~38 arcsec (3 kpc in projection) 
southeast of the nucleus of M 87, the central dominant galaxy of the Virgo Cluster. We report the detection of 
far-infrared (FIR) [C II] line emission at 158 /mi from the nebulae using observations made with the Herschel 
Photodetector Array Camera & Spectrometer (PACS). The infrared line emission is extended and cospatial 
with optical Ha+ [N II], far-ultraviolet C IV lines, and soft X-ray emission. The filamentary nebulae evidently 
contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~ 100 K 
to ~ 10 7 K. This material has most likely been uplifted by the Active Galactic Nucleus from the center of 
M 87. The thermal pressure of the 10 4 K phase appears to be significantly lower than that of the surrounding 
hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the 
filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance 
the pressure of the surrounding ICM is B ~ 30 - 70 fiG. The spectral properties of the soft X-ray emission 
from the filaments indicate that it is due to thermal plasma with kT ~ 0.5-1 keV, which is cooling by mixing 
with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. 
Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent 
with the apparent lack of X-ray emitting gas with kT < 0.5 keV. However, we show that the missing very soft 
X-ray emission could be absorbed by the cold gas in the filaments with an average absorption column density 
of A^h ~ 10 21 cm 2 , providing a natural explanation for the apparent temperature floor to the X-ray emission at 
kT ~ 0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles 
penetrating the cold gas following a shearing induced mixing process. An additional source of energy may, in 
principle, be provided by X-ray photoionization from cooling X-ray emitting plasma. The relatively small line 
ratio of [O I]/[C II] < 7.2 indicates a large optical depth in the FIR lines. The large optical depth in the FIR lines 
and the intrinsic absorption inferred from the X-ray and optical data imply significant reservoirs of cold atomic 
and molecular gas distributed in filaments with small volume filling fraction, but large area covering factor. 
Subject headings: galaxies: individual (M 87) - galaxies: clusters: intracluster medium - infrared: ISM 



1. INTRODUCTION 

Massive giant elliptical galaxies in the centers of clusters 
with central cooling times shorter than the Hubble time fre- 
quently display spectacular optical emission-l i ne nebulae (e.g. 
Johnstone etlillT98"7l JHeckman et aT1[T989l |Donahue et al.| 
1992| |Crawford et al.||1999| |McDonald et al.||2010) . Asso- 
ciated cold molecular gas has also been detected in many of 
these systems, from near-infrared (NIR) K-band H2 fine emis- 



et al. 2000 1 Hatch et al. 2005; 



sion (e.g.lJaffe & Bremer|1997l|Falcke et al.|1998||Donahue 

"i Jafi'e et al.||2005i J ohnstone 



et al.|2 007 ; Oonk et al. 2010) a nd CO observations ( e.g. Edge 
2T)jTTp :dge & Frayer 2003} |Salome & Combes|26^|Mcbon- 
ald et al.|2012| >. While the NIR spectra show H2 line ratios 
characteristic of collisionally excited 1000-2000 K molecu- 
lar gas, the CO emission traces the coldest (<50 K) molecu- 
lar gas phases. Optical Ha+[N II] line emission arises from a 
thin ionized skin of 10 4 K gas on underlying reservoirs of cold 
neutral and molecular gas. CO observations indicate th at the 
amounts of molecular gas in cluster cores can be large (Edge 



™ Ha+ 



2001 ), e.g. the cold gas mass in the center of the Perseus Clus- 
teris approaching 10 11 M Q ( jSalome et al.|2006||2008j|201 1| ). 
The observed level of star formation in these systems is, how- 
ever, typically relatively small (e.g. |Q'Dea et al.|2008] l. De- 
spite significant efforts to understand the origin and the excita- 
tion mechanism of these nebulae, their nature, energy source, 

and detailed physics remain poorly understood; 

The relatively nearby (d -16.7 Mpc, Blak eslee et al.|2009) > 
giant elliptical galaxy M 87 (NGC 4486), at the centre of 
the Virgo Cluster, harbors a well-known extended optical 
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N II] filament system ( [Ford & Butcher|[T979t |Sparks 
of 



1993 ). The proximity of this system allows us to study 



lysics of the emission-line nebulae in greater detail than 
is possible elsewhere. The Ha+[N II] nebulae in M 87 were 
found to spatially coincide with filamentary s oft X-ray emis- 
sion ( |Young et al.||2002| Sparks et al. 2004), which can be 
modeled " as ~0.5 keV plasma in collisional ionization equi- 
librium ( |Werner et al.|2010| ). The emission-line nebulae also 
spatially coincide with C IV line emission at far-ultraviolet 
wavelengths (FUV), which typically arises in gas at temper- 
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TABLE 1 

Summary of observations and of the properties of the FIR lines integrated over the full 5x5 spaxel (47"x47") field of view of Herschel PACS. 



Line 


Rest frame A 


Observation duration 


Line Flux 


Instrumental FWHM 


Observed FWHM 


Line shift 




(fim) 


(s) 


(10 -14 erg s- 1 cm -2 ) 


(kms" 1 ) 


(kms- 1 ) 


(kms- 1 ) 


Cn 


157.7 


3000 


5.13 ±0.5 


240 


361 ± 27 


-62 ± 11 


O I 


63.2 


3312 


<3fi.fj 


85 






Olb 


145.5 


4480 


< 2.9° 


255 







a 2cr upper limits 

ature ~10 5 K ( |Sparks et al |2009|[2012] l. C IV line emission 
has so far not been detected in extended filaments in any other 
central cluster galaxy. Wern er et al.| (|2010) showed that all 
of the bright Ha and UV filaments in M 8/ are found in the 
downstream region of a <3 Myr o ld shock front revea led by 
X-ray observations with Chandra (Million e t al.|2010| . This 
argues that the generation of Ha, UV, and soft X-ray emission 
in M 87 is coupled to shocks in the hot X-ray emitting plasma. 
Based on these observations Werner et al. (2010) proposed 
that shocks induce shearing around the cooler, denser gas fil- 
aments, which promotes mixing with the ambient hot intra- 
cluster medium (ICM) via instabilities. By enhancing the 
rate at which hot thermal particles come into contact with the 
colder gas phases, mixing can in principle supply the power 
and the ionizing particles needed to explain the NIR to FU V 
line emission ( [Ferland et al.|2008[ [20091 |Fabian et al.|2011t . 

The launch of the Herschel space observatory, with its 
unprecedented sensi tivity to far-infrared (FIR) line emission 
(Pilbr atFet al.|2010| ), has opened new opportunities to study 
the coldest gas phases in galaxies. Herschel has enabled the 
first detections of atomic cooling lines of [C II], [O I], and 
[N II] in the X-ra y bright cores of Abell 1068, Abell 2597 
(Edge eT al.||2010[ ) , and the Centaurus and Perseus clusters 
(Mitta fet al.|201 l[|2012| l. These lines are the dominant cool- 
ing lines for neutral interstellar gas and can be used as diag- 
nostics toinferterr^eratures, densities, and radiation fields 
(e.g. |Kaufman et al.|19 99). Because carbon is the fourth most 
abundant element in the Universe and has a low ionization 
potential, the 158 fim [C II] line is the mos t ubiquitous and 
best-studied (e.g. |Malhotra et al . II 997] 1 200 1 ) > . It is a tracer of 
gas with a temperature of ~ 100 K. 

To study the cold gas phases associated with the filamen- 
tary line emission nebulae in M 87 we observed the regions 
of brightest Ha+[N II] filaments, extending to the southeast 
of the nucleus of the galaxy, with the Herschel Photodetec- 
tor Array Camera & Spectrometer (PACS) at the wavelengths 
of [C II]/ll57 y um, [O I]/i63yum and [O Ib]^145//m. Here we 
report the results of these observations, as well as a reanal- 
ysis of deep Chandra X-ray, Hubble Space Telescope (HST) 
optical and UV data, and long slit spectra obtained with the 
William Herschel Telescope (WHT). Sect. |2]describes the ob- 
servations, data reduction, and analysis of the Herschel PACS, 
HST, Chandra, and WHT data. In Sect. [3] we summarize the 
results, and in Sect. [4] discuss the implications of our obser- 
vations for the nature and the energy sources of the filaments. 
Our main conclusions are sum marized in Sect. [3 We assume 
a distance to M 87 of 16.7 Mpc ( jBlakeslee et al.|2009[ ), which 
implies a linear scale of 8 1 pc arcsec" 1 . The redshift of M 87 
is z = 0.004360. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Far Infrared spectroscopy with Herschel PACS 

We observed the FIR cooling lines of [C U]Al51fim, 
[O I]A63/*m, and [O Ib]/H45/xm with the PACS integral-field 
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Fig. 1. — The FIR [C II]/1157 fim line obtained from the central spaxel 
(9.4 X 9.4 arcsec ) of the Herschel PACS rebinned data cube. The mean 
line centroid is blue-shifted with respect to M 87 (at z = 0.004360) by 
v = - 123 ± 5 km s . The blue line indicates the detector noise. 

spectrometer (Poglitsch et al. 2010) on the Herschel space ob- 
servatory (Pilbr att et al.|2010| . The observations, with a dura- 
tion of 10880 seconds, were performed on January 4th, 2012 
(Obs. ID: 1342236278). Table[T]gives a summary of the ob- 
servations. Columns list the observed lines, their rest frame 
wavelengths, observation durations, the observed line fluxes 
and the 2cr upper limits, the spectral resolution full-width-at- 
half-maximum (FWHM) of the instrument at the given wave- 
length, the observed FWHM, and the observed line shift with 
respect to the systemic velocity of M 87 (v = 1307 km s ). 

The observations were taken in line spectroscopy mode 
with chopping-nodding to remove the telescope background, 
sky background and dark current. A chopper throw of 6' was 
used. All three line observations were taken in pointed mode 
centered on (a, <5)=(12:30:51.5, +12:23:07) (J2000). 

The observations were reduced using the HIPE software 
version 8.2.0, using the PACS ChopNodLineScan pipeline 
script for pointed observations. This script processes the data 
from level (raw channel data) to level 2 (flux calibrated spec- 
tral cubes). 

During the final stage of reduction the data were spectrally 
and spatially rebinned into a 5x5x/l cube. In the following 
we will refer to this cube as the re-binned cube. Each spatial 
pixel, termed spaxel, in this cube has a size of 9.4"x9.4". The 
cube thus provides us with a field of view (FoV) of 47"x47". 

For the wavelength re-gridding we set the parameters over- 
sample and upsample equal to 2 and 1 respectively. This 
means that one spectral bin corresponds to the native reso- 
lution of the PACS instrument (see the PACS Data Reduction 



Filamentary cold gas in M 87 



3 




187.722 187.720 187.718 187.716 1B7.714 187.712 187.710 
Right ascension 




[CM] velocity 



1 87.722 187.720 1 87.71 8 1 B7.716 187.71 4 1 87.71 2 187.71 
Right a 




187.722 187.720 187.718 187.716 187.714 187.712 187.710 
Right ascension 



-150 -100 



1 00 1 50 



: -' -1 kpc ■ ''■ _ 


/ 




■^'■Nprthern filament J 










Southern filament 












* . ' ' '. is 1 




Ha from HST WFPC2 








FUV from HST ACS/SBC 



187.722 187.720 187.718 187.716 187.714 187.712 T&rVT+fl-' 
Right ascension 



187.722 187.720 187.71B 1B7.716 187.714 187.712 187.710 
Right ascension 




187.722 187.720 187.718 187.716 187.714 187.712 187.710 
Right a 



1 kpc 



t 



0.5 keV gas 



187.722 187.720 187.718 1B7.716 187.714 187.712 187.710 
Right ascension 




187.722 187.720 1B7.718 187.716 187.714 187.712 187.710 
Right ascension 



20 30 40 50 60 70 




200 300 400 



1B7.722 187.720 187.718 187.716 187.714 187.712 187.710 
Right ascension 



400 500 600 700 800 



Fig. 2. — FIR [C II] maps (top panels) and Ha-+[N II], FUV, and X-ray images of the filamentary multiphase gas southeast of the nucleus of M 87. All images 
show the same 47 X 47 arcsec region of the sky. The projected distance of the bright filaments from the nucleus is 38 arcsec (3 kpc). Top left panel: Map of 
the integrated [C II] line flux in units of 10~' 4 erg s -1 cm~ 2 per 6"x6"spaxel obtained with Herschel PACS. Top central panel: The velocity distribution of the 
[C II] emitting gas, in units of km s -1 , relative to the systemic velocity of M 87, v = 1307 km s . Top right panel: Map of the velocity dispersion, cr, of the [C II] 
emitting gas. Middle row left panel: Ha+[N II] image obtained with HST WFPC2. The contours of the se filaments are over-plotted on the Herschel PACS maps 
in the top row. On this HST WFPC2 image, contours of 6 cm radio emission from |Hines eTaT]fT989| are over-plotted in red. Middle row central panel: FUV 
image showing C IV line emission obtained with HST ACS/SBC. Middle row right panel: Chandra sort X-ray image, extracted in the 0.7-0.9 keV band. Bottom 



i right panel 

panels: Spatial distribution of the emission measure, Y = J niin^dV, (in 10 58 cirT 3 
(bottom right) plasma detected at 99.7 per cent confidence. 



arcsec ) of the 0.5 keV (bottom left), 1.0 keV (bottom central), 2.0 keV 
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GuideQfor further information). Larger values for both these 
two parameters were investigated and did not change the re- 
sults. 

For the [C II] line we projected the resinned cube onto 
the sky using the specProject task in HIPE and the hrebin 
task in IDL. In the following we will refer to this cube as the 
projected cube. This projected cube has a pixel size of 6"x6" 
and the same FoV as the resinned cube. 

2.2. HST Ha+[N II] and FUV data 

To study the detailed morphology of the filamentary 
emission-line nebulae, we have also analyzed Ha+[N II] im- 
ages taken with the Wide Field Planetary Camera 2 (WFPC2) 
through the F658N filter for 2700 seconds (proposal ID: 
5 122). The central wavelength of the filter is at 6591 A and its 
bandwidth is 29 A. It transmits the Ha line at A = 6563A, with 
transmission T = 0.80, and two [N II] lines at A = 6584A, 
with T = 0.20 and at A =6548A, with T = 0.76, at the sys- 
temic velocity of M 87 (v = 1307 km s" 1 ). We subtracted the 
emission of the underlying stellar population of the galaxy 
using a red continuum image taken with WFPC2 through the 
filter F8 14W (proposal ID: 5941). 

A FUV image was obtained using the Advanced Camera 
for Surveys Solar Blind Channel (ACS/SBC) (proposal ID: 
1 1861) through the filter F150LP (1 630 seconds) whic h cov- 
ers the C IV line at A = 1549 A (see|Sparks et al.|2009|>. 



2.3. Chandra X-ray data 

Extensive Chandra X-ray observations of M 87 were made 
between July 2002 and November 2005 using the Advanced 
CCD Imaging Spectrometer (ACIS). The total net exposure 
time after cleaning is 574 ks. The data reduction is described 



in [Mill ion et al. ( |2010| l. We extracted six background sub- 
tracted, flat-fielded, narrow band images between 0.3 keV and 
2.0 keV with a 0.492 x 0.492 arcsec 2 pixel scale (raw detector 
pixels). To measure the detailed properties of the soft X-ray 
emission, we also extracted spectra fro m regions dete rmined 
using the Contour Binni ng algorithm (Sanders 2006), as de- 
scribed in Werner et al. (2010). Spectral modeling has been 
performed with the SPEX package ([Kaastra et al. 1996) in the 
0.5-7.0 keV band. 

2.4. Long slit optical spectra 

Long slit spectroscopy was performed using the Interme- 
diate dispersion Spectrograph and Imaging System (ISIS) at 
the 4.2 m William Herschel Telescope on the island of La 
Palma on February 10, 2011. The total on-source integra- 
tion time was 3500 s, split into 7 exposures, for the blue and 
red arm simultaneously. The data were reduced using stan- 
dard IRAF (Image Reduction and Analysis Facility) proce- 
dures. The CCD images were bias subtracted, flat fielded, 
wavelength calibrated, flux calibrated, and merged. Spectra 
were shifted into the heliocentric velocity frame. Cosmic rays 
were removed with the CRNEBULA subroutine, designed for 
diffuse objects. Blue and red arm images were treated sepa- 
rately during the whole analysis. 

3. RESULTS 
3.1. Cold and warm gas phases 

1 http://herschel.esac.esa.int/twiki/pub/Public/PacsCalibrationWeb/ 
PDRG_Spec_Mayl2.pdf 



We detect the filamentary emission-line nebulae across all 
wavebands studied, from the FIR to the soft X-rays. Using 
Herschel PACS, we detect [C U]Al51pm line emission (see 
Fig.[T} and determine 2<x upper limits for the [O I]/163^/m and 
[O Ib]/I145^m cooling lines. The properties of the FIR lines, 
spatially integrated over the 5x5 spaxel (47"x47") Herschel 
PACS field of view, are summarized in TablefT] The 2cr upper 
limits on the fluxes of the [O \\A63pm and[0 Ib]/ll45/im 
lines given in Table [T] were determined assuming that their 
velocity widths are equal to the FWHM of the [C II] line. 

The [C II]/H57yum line emission is extended and spatially 
coincident with the Ha+fN II] filaments (see the top left 
panel of Fig.|2]i. We see two prominent filamentary structures, 
which we win call the 'southern filament' and the 'northern 
filament'. The brightest southern filamentary structure falls 
on the central spaxel of the PACS spectrometer. The [C II] 
flux detected from the central spaxel (9.4 x 9.4 arcsec 2 ) in the 
rebinned cube is (8.5 + 0.7) x 10~ 15 erg s" 1 cirT 2 . The top 
central panel of Fig. [2] shows a map of the velocity distribu- 
tion of the [C II] line emitting gas with respect to the velocity 
of M 87. This image shows an interesting dichotomy, with 
the northern filaments receding at about +140 km s" 1 and 
the bright southern filament (most of the emission of which 
falls on the central 9.4 x 9.4 arcsec 2 spaxel of the detector) 
moving in the opposite direction along our line of sight at 
-123 ±5 km s" 1 , with respect to M 87. These velocities show 
good agreement with the radial velocities of the correspond- 
ing optical filament s measured using long-slit spectroscopy 
( Spark s et al.|Jl99 3"). The top right panel of Fig. [2] shows a 
map of the velocity dispersion of the [C II] line emitting ma- 
terial. While the northern, receding filaments show a signifi- 
cant velocity dispersion of cr = FWHM/2.355 ~ 124 km s" 1 

(FWHM= V377 2 - 240 2 = 291 km s _1 , corrected for the 
instrumental resolution FWHM of 240 km s" 1 ), the velocity 
dispersion measured in the southern filaments, cr ~ 55 km s" 1 

(FWHM= V273 2 - 240 2 = 130 km s" 1 ), is just marginally 
above the instrumental resolution. 

The middle row left panel of Fig. [2] shows the Hq-+[N II] 
image obtained with HST W FPC2. Contours of 6 cm radio 
emission from Hines et al. ( 1989) are over-plotted, showing 
that in projection the side of the southern filament closest 
to the nucleus overlaps with the radio lobes. The narrow- 
est resolved filaments of Har+[N II] emission have a diam- 
eter of only 0.4 arcsec, which corresponds to 32 pc. The total 
count rate measured in our narrow band image from the bright 
southern filamentary structure within a 15.7"x5.7" box- 
like region centered at (a,J)=(12:30:51.48,+12:23:07.33) is 
6.94 counts s" 1 . Assuming an [N II]/16548/Hq' flux ratio of 
0.81 and [N n],i6584/Har=2.45 ( |Ford & Butcherl[T979l > and 
taking into account the filter throughput at the wavelengths of 
the lines (see Sect. I2.2[ ), the Ha flux of the southern filament 
is faa = 9.0 X lO -1 ^ erg s" 1 cm 2 . 



Sparks et al. ( 2009b reported that the Ha+[N II] emission- 
line nebulae are co-spatial with FUV emission (see the mid- 
dle row central panel of Fig. [2ji. Using subsequent spectro- 
scopic measurements with the HST Cosmic Origin Spectro- 



graph, Sparks et al. ( 2012} showed that the FUV emission in 
the northern filaments is due to C IV line emission. We find 
that the C IV flux, measured within the same aperture as the 
Ha+[N II] flux above, is few - 1-2 x 10" 14 erg s" 1 cm" 2 . 

The optical spectral properties of the southern filament 
show an intriguing rapid decrease in the [O III]/Hj6 ratio with 
increasing distance from the nucleus, dropping monotonically 
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Fig. 3. — Left panel: Optical diagnostic diagram (a BPT diagram, Baldwin et al. 1981 1 showing the [O III]/H/? against the [N Ilj/Hcr flux observed in the southern 
filament at different radii from the nucleus. The brightening shades of gray indicate increasing projected distance from the nucleus in arcsec. While the side of 
the filam ent closest to the nucl eus has spectra similar to Seyfert l ike AGN, at larger radii the gas has a LINER-like emission spectrum. The solid and dashed lines 
are from Kewley et al. 1 2006 1 and Osterbrock & Ferland (2006 1, respectively. Right panel: The measured [S II]/16716/[S II]/16731 line ratios, a good probe of 
gas density, as a function of radius. T he full, dashed, and dotted lines indicate ratios corresponding to electron densities of n t = 30 cirT 3 , n t = 100 cm' 3 , and 
n e = 500 cm , respectively jOsterbrock & Ferland|2006] . 



from ~2.5 at a projected distance of 28 arcsec to ~0.7 at 43 
arcsec (see Fig. [3]). The diagram in the left panel of Fig. J3] 
showing the [O III]/Hj6 over the [N ll]/Ua ratios (a BPT dia- 
gram, |Baldwin et al.|198"T| l, indicates that the side of the fila- 
ment closest to the nucleus, which also overlaps with the radio 
lobes (see middle row left panel of Fig. [2]), has spectra similar 
to Seyfert like active galactic nuclei (AGN). The [O III]/Hy6 
line ratios in this part of the filament are larger than in the nu- 
cleus of M 87, where we measure a ratio of ~ 1.75. At larger 
radii, the gas has a LINEB^J-like emission spectrum, which 
is more similar to the emission line spectra typically seen in 
extended filaments around brightest cluster galaxies in cluster 
cores. 

The right hand panel of Fig. [3] shows the measured 
[S II]/16716/[S II]/16731 line ratio, a good probe of gas den- 
sity, as a function of radius. The full, dashed, and dotted 
lines indicate ratios corresponding to electron densities of 
n e = 30 cm' 3 , » P _ 100 cm' 3 , and n e = 500 cm' 3 , respec- 
tively (Osterbrock & Ferland|2006). The value n e - 30 cm' 3 



is very close to the low density limit, below which the [S II] 
line ratios do not provide useful constraints on gas density. 
Along most of the filament, the line ratios are close to this 
limit indicating that the electron density of the 10 4 K phase 
is n e < 30 cm ' 3 . At radii between 29-30.5 arcsec, the mea- 
sured line ratios are outside of the low density limit. The data 
point closest to the nucleus, at 28 arcsec, indicates a possible 
increase in density, but with large error bars. 

3.2. The soft X-ray emission 
Young et al.] d2002) and |Sparks et ah] (|2004|l found good 



middle row right panel of Fig. E|. The spatial correlation be- 
tween the 0.7-0.9 keV X-rays and the Har+[N II] line emis- 
sion is remarkable: everywhere, where we see Ha+[N II] 
emitting gas, we also see an excess of soft X-rays. On the 
other hand, in the 0.5-0.7 keV and 0.9-1.2 keV bands the 
filaments do not stand out clearly against the surrounding 
emission, indicating that most of their X-ray emission comes 
from the Fe XVII and Fe XVIII lines in the 0.7-0.9 keV 
band. The excess X-ray flux of the filaments in this band 
is 1.1 x 10' 14 ergs s' 1 cm' 2 and their average X-ray surface 
brightness is 1.2 x 10' 16 ergs s' 1 cm' 2 arcsec' 2 . 

Assuming that the soft X-ray emission is of a thermal ori- 
gin, we follow the analysis of Werner et al. (2010! and fit 
to each spatial bin a model consisting of collisionally ion- 
ized equilibrium plasmas at three fixed temperatures (0.5 keV, 
1.0 keV, and 2.0 keV), with variable normalizations and com- 
mon metallicity. The bottom panels of Fig.|2]show the spatial 
distributions of the emission measures of the individual tem- 
perature components. The 0.5 keV component closely fol- 
lows the distribution of the cold gas phases. This compo- 
nent is required to fit the Fe XVII and Fe XVIII lines in the 
0.7-0.9 keV band, which are only present in regions where 
Ha+[N II] emission has also been detected. The spatial cor- 
relation between Ha+[N II] and the softest X-ray emitting 
component is, however, not perfect: the ratio of Ha+[N II] 
flux to the ~0.5 keV emission component is ~ 5 times larger 
in the southern than in the northern filament. The 1 keV com- 
ponent is enhanced as well in the regions where cold gas is 
present, but it is spatially more extended and is also seen in 
many regions where cold gas is not detected (see|Werner et al. 



spatial correlation between Ha+|N II] fine emission and soft |2010| l. The spatial distribution of the 2 keV plasma does not 
X-ray emission. Here, we have analyzed almost four times as 
much Chandra data, which provides excellent photon statis- 
tics in the regions of interest. We extracted six narrow band 
images between 0.3 keV and 2.0 keV and found that the fil- 
aments stand out clearly in the 0.7-0.9 keV band (see the 



Low-ionization nuclear emission-line region ( Heckman 1980 



correlate with the Ha+[N II] emission. 

Because the soft X-ray line emission observed in the fila- 
ments may be due to cooling gas, we fitted the spectra with 
a model consisting of a single-temperature plasma in colli- 
sional ionization equilibrium plus an isobaric cooling flow, 
which models cooling between two temperatures at a certain 
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metallicity and mass deposition rate. The upper temperature 
and the metallicity of the cooling flow component were tied 
to the temperature and metal abundance of the single- temper- 
ature plasma. Both components were absorbed by a Galactic 
absorption column density Nu = 1.9 X 10 2l) crrT 2 (IKalberla 
|et al.|2005j >. Fixing the lower temperature cut-off of the cool- 
ing flow model to kT\ 0VJ = 0.5 keV, we obtained a best fit up- 
per temperature of kT up = 1 .90 ± 0.04 keV, a mass deposition 
rate of M = (2.46 + 0.13) x 10~ 2 M e yr" 1 , and a met allicity of 
Z - 1.46+ 0.09 Solar (relative to the Solar values by Grevesse 
& Sauval|1998| . Even though the spectra indicate that above 
kT ~ 0.5 keV the gas may be cooling, the spectral signatures 
of gas cooling below 0.5 keV - soft X-ray emission, includ- 
ing the O VII line - are missing (see also the high resolution 



spectra in [Werner et al. 12006 p010| >. In order to test whether 
the missing soft X-ray flux might have been absorbed by the 
cold gas in the filaments, we fixed the mass deposition rate in 
the model to the best fit value obtained with a lower temper- 
ature cut-off set to 0.5 keV, extended this cut-off down to a 
1000 times lower temperature of 0.5 eV (the lowest value al- 
lowed by the model), and included an intrinsic absorber in our 
spectral model. This model provides a good fit with intrinsic 
absorption of N n = (1-24 ± 0.12) x 10 21 cm ~ 2 and absorbed 
bolometric flux of 3.6 x 10~ 14 erg s" 1 cm -2 . 

4. DISCUSSION 

4.1. Magnetized filaments of multi-phase material 

In the filamentary emission-line nebulae to the southeast 
of the nucleus of M 87, we detect gas spanning a tempera- 
ture range of at least 5 orders of magnitude, from ~ 100 K 
to ~ 10° K. The [C II], Ha+[N II], and C IV emitting gas 
phases are consistent with being co-spatial, forming a multi- 
phase medium. Soft X-ray emission in the 0.7-0.9 keV band 
is also always present in regions where Ha+ [N II] emission 
is seen. Its spectral shape is consistent with line emission of 
~0.5 keV thermal plasma, previously detected in the XMM- 
Newton high resolution reflection grating spectra integrated 



over the central region of M 87 (|Werner 



etaD : 
s with 



2006 



2010). 



We see no evidence of X-ray emitting gas with temperature 
kT < 0.5 keV. 

The [C II] 158 pm line is about 1500 times stronger than 
the CO (1— »0) rotational line in normal galaxies and Galac- 
tic molecular clouds, and 6300 times more intense in star- 



burst galaxi e s and Galactic star forming regions (Crawford 
it~aT1[T985| |Stacey et al^|T99T) . Given a [C II] flux of 
8.5 X 10~ ls erg s" 1 cm~ z observed from the brightest region 
of the southern filament, the expected CO (1— >0) flux is 1.3 
5.3 x lO' 18 



erg s 



cm 



2xlQ-' 7 erg s ' T cm' 2 for CO (1 



2 consistent with the upper limit of 
0) at 2.6 mm observed in the 



same region (Salome & Combes 2008). Using the standard 
CO luminosity to H2 conversion factors the inferred range of 
CO luminosities corresponds to an H2 mass ofO.4-2xlO 6 M . 

Assuming that all the observed gas phases are in thermal 
pressure equilibrium with the ambient ICM (assuming p = 
0.1 1 - 0.22 keV cirT 3 , for a distance range of r = 3 - 6 kpc 
from the nucleus depend ing on the pos ition of the filament 
along our line of sight, |Churazov et al.||2008| ), the ~ 100 K 
[C II] emitting gas forms a network of narrow filaments with 
densities of the order of (1.3 -2.6)x 10 4 cm 3 (higher than the 
critical density of 3 x 10 3 cm" 3 above which the gas is in local 
thermal equilibrium with the level populations determined by 
collisions) with volume filling fraction of fy ~ (1 - 2) x 10~ 4 
in the bright southern filament. Although the presence of O I 



6300 A lines in the optical spectra (see also Ford & Butcher 
1979 ) indicates the existence of large partially ionized regions 
in the filaments, the presence of [S II] lines most probably 
indicates a fully ionized ~ 10 4 K phase surrounding the colder 
gas. The density of this H<?+[N II] emitting gas, assuming 
thermal pressure equilibrium with the ambient ICM, should 
be - (1.3- 2.6) x 10 2 cirT 3 . 

The [S n]/l6716/[S HJA6731 line ratios, however, indi- 
cate much lower particle densities in the 10 4 K phase (n ~ 
2n e < 60 crrT 3 ) implying the presence of additional signifi- 
cant non-thermal pressure components in the filaments, such 
as turbulence and magnetic fields. Assuming isotropic micro- 
turbulence with a characteristic velocity Vturb m a medium 
with a sound speed c s , the turbulent pressure support will 
be p ta rb = l/2y/?therm M 2 , where y is the adiabatic index 
assumed to be 5/3 and M = v tur b/c s is the Mach number. 
However, assuming characteristic turbulent velocities of the 
order of the sound speed in the 10 4 K phase (c s = 15 km s" 1 ), 
the sum of thermal and turbulent pressure, /? tur b + ptherm = 
0.09 keV cirT 3 (1.4 x 10~ 10 dyne cm -2 ), is still smaller than 
the surrounding ICM pressure. If the turbulence in the Ha 
emitting phase is subsonic, then the additional pressure will 
be provided by magnetic fields. A magnetic pressure of 
Pmag = B 2 /8tt ~ (0.3 - 2.1) x 10" 10 dyne cnr 2 , needed to 
keep the filaments in pressure equilibrium with the surround- 
ing ICM, requires magnetic fields of B — 28 - 73 pG. This 
is similar to the values inferred using arguments based on 
the integrity of Hg +[N II] filaments in the Perseus Cluster 
(Fabian et al. 2008 ). Radio observations of the Faraday rota- 
tion measure in a number of coolin g cores revealed magnetic 
field strengths of 10-25 pG (Taylor et al.|200T] [2007) [Allen 



et al.|200T Feretti et al.|1999 1. Significant imbalance in ther- 
mal pressure of the gaseous filaments in cool-core clusters 
has previously also been inferred between the warm molec- 
ular hydrogen gas seen in NIR (nT - 10 8 ~ 9 cm ~ 3 K) and the 
ionized gas phases seen in the optical (nT ~ 10 5 cm" 3 K), 
emphasizing the need for dyna mic models for the filaments 
Paffe et al.|200T]|Oonk et al.|2010p . 

The smallest resolved diameter of the Ha+[N II] filaments 
is 32 pc, about half of the 70 pc diameter observed in the 
more distant Perseus Cluster (Fabi an et al.||2008). However, 
even the narrowest resolved filaments are likely to consist of 
many smaller strands. They are surrounded by hotter 10 5 K 
C IV emitting gas. The space between these filaments is filled 
by X-ray emitting gas. If the soft X-ray line emission is due 
to ~0.5 keV gas, then based on our three-temperature fits (see 
Sect. [3} and assuming thermal pressure equilibrium with the 
ambient ICM, its mass in the southern filament is ~ 10 5 M Q 
and its emitting volume is ~0.016 kpc 3 . 

4.2. The energy source of the filaments 

Spar ks et al.| ([2009 ) show that no main sequence stars of 
early spectral type are present in the filaments and therefore 
neither the FUV nor the Ha+[N II] emission can be due to 
photoionization by hot UV emitting stars. The [O III]/Hy0 
emission line ratios are high on the side of the southern fil- 
ament closest to the nucleus, which also overlaps with the 
radio lobes (see middle row left panel of Fig. [2]), and then 
decrease rapidly with increasing distance from the nucleus. 
The fact that the line ratios in the filament reach values which 
are higher than those measured in the nucleus argues against 
photoionization by the AGN. It seems more likely that the in- 
creased [O III]/Hy6 ratios are due to internal shocks produced 
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as the radio lobes push against the cold gas. These shocks 
seem to affect only a part of the filament and the emitted Ha 
flux does not decrease with the decreasing [O III]/Hy8 emis- 
sion line ratios. Therefore, they can be ruled out as energy 
sourc e powering the emission of these nebulae. 

Ferl and et al.| ( |2009] l showed that the broad band (optical to 
mm) emission-line spectra of nebulae around central galaxies 
of cooling core clusters most likely originate in gas exposed 
to ionizing particles. Such ionizing particles can be either rel- 
ativistic cosmic rays or hot particles penetrating into the cold 
gas from the surrounding ICM. Saturated conduction from the 
hot into the cold phase of the southern filament would produce 
a heat flux of 5(ppc s ~ 0.1 erg s _1 cirT 2 , where <p ~ 1 accounts 
for the uncertain physics, p = 0.22 keV cm ~ 3 is the gas pres- 
sure a nd c s =510 km s~' is the sound speed in the hot 1 keV 
phase ( |Cowie & McKee][T977l Fabian et al.pOTT] ). Follow- 
ing the argument of [Fabian et al.| ( |2011 1, given a typical Ha 



surface brightness of ~ 2 x 10~ lb erg s _1 cm" 2 arcsec -2 , cor- 
responding to an emitted flux of 1.1 x 10~ 4 erg s _1 cirT 2 and 
assuming that the total broad band emitted flux is 20 times 
larger 2.2 x 10~ 3 erg s ~' cirT 2 , the required efficiency with 
which impinging thermal particles penetrate and excite the 
cold gas is of the order of / = 0.022 (/ = 0.044 for a fil- 
ament at r — 6 kpc from the nucleus), which is reasonable. 
However, these particles must somehow overcome the obsta- 
cle presented by the magnetic fields believed to be at the in- 
te rface of the cold and hot gas phases. 

| Werner et aT] ( |2010| ) show that all of the bright Ha+[N II] 
and UV filaments are seen in the downstrea m region of a 
< 3 M yr old AGN induced M > 1 .2 shock front (Milli on etaE] 
20 IP} that seems to have just passed the filaments. Based on 



these observations they propose that shocks induce shearing 
around the cold, filamentary, dense gas, thereby promoting 
mixing of t he cold gas with the ambient hot ICM via instabil- 
ities (e.g. Fried man et al.|2012[ ). Instabilities and turbulence 
may also be driven by shearing motions between the hotter 
X-ray gas and the filaments, as they move through the ambi- 
ent ICM, both as they are being uplifted and as they fall back. 
This type of more continuous shear induced mixing could per- 
haps be in the long term more important for energizing the fil- 
aments than shocks. Fabi an et al.| (f2"0l 1) propose that the sub- 
mm through UV emission seen in NGC 1275 in the core of the 
Perseus Cluster is due to the hot ICM efficie ntly penetrating 
the cold gas through reconnection diffusion (|Lazarian et al. 



|2010 201 1 1. The same process may also be operating in the 



filaments of M 87 with the fast reconnection induced by shear- 
ing instabilities and turbulence. Shearing and magnetic recon- 
nection are, however, not essential to bring the cold gas into 
contact with the hot ICM. The support from magnetic fields 
tied to the surrounding ionized phase may become unstable if 
the layer of dense cold neutral gas is deep enough, and there- 
fore it may fall through into direct contact with the hot gas. 
The hot ICM particles penetrating the cold gas then produce 
secondary electrons th at excite the ob s erved FIR through UV 
emission as shown by Ferland et al.| ( |2009[ ). Mixing of the 
cold and hot gas also naturally account for the presence of 
the ~ 10 5 K intermediate-temperature gas phase (geometric 
mean of the hot and cold gas temp eratures) predicted by the 
mixing layers calculated by Begelman & Fabian ( 1990| l, al- 
though mixing layers have most probabl y a m ore complicated 
temperature structure (Esquivel et al. 2006). As the X-ray 
emitting gas cools through mixing, the mass of the southern 
filament may grow by as much as 4 X 10~ 2 M Q yr _1 . 



The particle heating model of Ferland et al. (2009), calcu- 
lated considering optically thin emission from a cell of gas, 
predicts [O I]/[C II] ~ 21 and is inconsistent with the 2<x 
upper limit of 7.2 observed in M 87. The ratios observed in 
the Perseus and Centaurus clusters, and in Abell 1068 and 
2597 are even lower. This discrepancy has been discussed in 
a recent study of NGC 1275 by |Mittal etaL] ( |20 1 2\ , who con- 
clude that the line ratio suggests that the lines are optically 
thick, implying a large reservoir of cold atomic gas, which 
was not accounted for in previous inventories of the filament 
mass. The same can be concluded for M 87, where the large 
optical depth of the FIR lines (Canning et al. in prep.) and the 
intrinsic absorption inferred from the X-ray data (see below) 
imply significant reservoirs of cold atomic and molecular gas. 

4.3. X-ray line emission from cold gas due to charge 
exchange and inner shell ionization ? 

As ions from the hot ICM penetrate into the filaments, they 
pick up one or more electrons from the cold gas. The im- 
pinging ions thus get into excited states and their deexcitation 
may produce X-ray radiation. Astrophysical X-ray emission 
due to this 'charge excha nge' process was first discovere d in 
the comet Hyakutake ( |Lisse et"aL||1996| |Cravens||l997| ) and 
since then it has been identified as the mechanism res ponsi- 
ble for X-r ays from planets (e.g. Gladst one et al.|2002[|Den-| 
2002|), the geocorona (Wargeli n et al.||2 004; Fujimoto 
2007'), the heliosph ere (Robertson et al. 120011), and in- 



nerl 



et a 



terstellar medium in particular at interfaces between hot gas 
and cool clouds (see the review by Dennerl 2010 and ref- 
erences therein). Recently, charge exchange has also been 
found to possibly contribute significantly to the soft X-ray 



line emission from star-forming galaxies (Liu et al. 2011 



|2012| ). |Fabian et al.| ( |201 1) conclude that in the Perseus Clus- 
ter charge exchange may account for a few per cent of the 
soft X-ray emission from the filaments. They also speculate 
that if the filaments are composed of many small strands, the 
same ion might pass from the hot into the cold and back to the 
hot phase repeatedly, getting collisionally ionized while pass- 
ing through the hot ICM, thus enhancing the charge exchange 
emission to the observed level. 

After hydrogen and helium, the most abundant ions pene- 
trating the cold gas are O s+ , C 6+ , N 7+ , Ne 10+ , Mg 12+ , Mg I1+ 
(with concentrations 8.4 x 10~ 4 , 3.6 x 10~ 4 , 1.1 x 10~ 4 , 
1 . 1 x 10~ 4 , 2.5 x 10~ 5 , 1 .07 x 10~ 5 relative to hydrogen by num- 
ber). Charge transfer onto these ions is expected to produce 
mainly hydrogen and helium like line emission in the X-ray 
band. In the 1 keV plasma surrounding the filaments, 23% of 
Fe is in the form of Fe 20+ , another 23% in Fe 21+ , 18% inFe 19+ , 
and 17% in Fe 22+ , with relative concentrations smaller than 
1.1 x 10~ 5 . The charge exchange line fluxes produced by Fe 
will therefore be significantly lower, with most of the line flux 
in Fe XIX to Fe XXII lines, in the spectral band above 0.9 keV. 
Moreover, laboratory measurements of X-ray spectral signa- 
tures ofcharge exchange in L- shell Fe ions by Beiersdorfer 
|et al.| ( |2008| ) found that the n > 4 — » n — 2 transitions are 
strongly enhanced relative to n = 3 — > n = 2, shifting the 
line emission to higher energies. Because charge transfer pro- 
cesses and the corresponding X-ray emission are expected to 
occur inside of the dense cold clouds, absorption may signifi- 
cantly decrease the observed fluxes of the lowest energy X-ray 
lines such as C VI, N VII, O VII, and O VIII. The Ne IX, Ne X, 
Mg XI, and Mg XII lines in the 0.9-1.9 keV band are, how- 
ever, not expected to be affected by absorption strongly and 
their enhancement might indicate charge exchange. Given 



8 



Werner et al. 



that none of these lines are enhanced in the filaments, and 
no strong charge exchange emission is expected in the 0.7- 
0.9 keV band, we conclude that the observed soft X-ray emis- 
sion is not due to charge transfer onto the impinging ions. 

The ions and electrons penetrating from the surrounding hot 
ICM into the cold filaments will also ionize the inner electron 
shells of the atoms and molecules in the cold clouds. This 
inner shell ionization may in principle also contribute to the 
X-ray emission of the filaments, but given the energetics, it is 
not expected to be significantly stronger than the X-ray emis- 
sion due to charge exchange. We would primarily expect to 
see X-ray lines from few times ionized C, N, O, Ne, and Mg, 
with the Ne and Mg lines in the band above 0.6 keV. Detailed 
predictions of this process are beyond the scope of this paper. 

Given that the spectral shape of the filaments is consistent 
with line emission of ~0.5 keV thermal plasma previously re- 
solve d within the central region of M 87 with XMM-Newton 
RGS ( ] Werner et al. 2006, 2010 ) and that charge exchange as a 
significant source of X-ray emission from the filaments can be 
ruled out, we conclude that the observed soft X-ray emission 
is most likely due to thermal plasma. We note, however, that 
although charge exchange is not responsible for the observed 
excess X-ray emission in the 0.7 - 0.9 keV band, the level 
which is required from the particle heating model is observa- 
tionally acceptable. 

4.4. AGN uplift induced cooling instabilities or mixing? 

The observed velocity distribution of the [C II] (see top cen- 
tral p anel of Fig. |2]l and optical emission lines (Spar kTet al.| 
1993) indicates that the northern filaments are currently being 
uplifted by the c ounter-jet (which i s oriented at < 19° from 
our line-of-sight, Biretta et al. 1999) moving away from M 87 
with a line-of-sight velocity of vlos = +140 km s _1 . The 
southern filaments, in contrast, are falling back towards M 87 
with v L os = -123 km s _1 . Assuming uplift at the Keplerian 
velocity of ~ 400 km s _1 , the southern filament would take 
8 x 10 6 yr to get to its current position. Given that this filament 
is most likely already falling back, its age is > 10 7 yr. Uplift 
of low entropy, metal enriched X-ray emitting gas by buoy- 
antly rising relativistic plasma from the AGN jets of M 87 has 
been previously inferred based on radio and X-ray data (|Chu-| 
razov et al.|2001HForman et al.|2005||Simionescu et al.|2008[ 



Werner et al. 2010[ ). Optical emission line nebulae also appear 
to have been drawn out by rising radio bubbles in other nearby 



systems (e.g. the Perseus a nd Centaurus clusters, Hatch et al. 
2006; CanningetamoTT) 



The initially uplifted material might have already been a 
multi-phase mixture of low entropy X-ray emitting plasma 



and dusty c old gas, as is observed in the core of M 87 ( |Sparks 
et al.|19 93 ). As the relatively dense, low entropy X-ray emit 



ting plasma was removed from the direct vicinity of the AGN 
jets it might have cooled radiatively, further contributing to 
the co ld gas mass of the uplifted material (see Reva z et al.| 
2008 ). The multi-phase X-ray emitting gas, with the 0.5 keV 



phase co-spatial with the Ho-+[N II] filaments and surrounded 
by a spatially more extended 1 keV plasma, indicates that such 
cooling may be taking place. Cooling instabilities are pre- 
dicted to take place in cooling core clusters within radii where 
the ratio of the cooling time to free fall time is f C ooi /fa $ 10 
( jSharma et al.|2012[|McCourt et al.|2012[|Gas pari et al .l2012| 
this condition is fulfilled at the radii where the M 87 fila- 
ments are seen). AGN induced uplift of multiphase gas may 
facilitate the onset of such local cooling instabilities. Some 
Hq'+[N II] and soft X-ray filaments lie at the edges of cavi- 



ties filled by radio emitting plasma, indicating that they might 
be sheets of gas seen in projection that have been swept up 
and compressed by the expanding and rising cavities. Even 
though such compression would first heat the gas, eventually 
it will accelerate its cooling. This picture of radiatively cool- 
ing X-ray gas is, however, complicated by the observed tem- 
perature floor at 0.5 keV. Several explanations were proposed 
to address the missing soft X-ray emission, including mixing 
(Fabian et al. 2 002|), anomalous O/ Fe abundance ratios, and 
absorption (Sanders & Fabia ripOl 1) . 

The mass deposition rate, inferred from the temperature 
distribution of the X-ray gas in the southern filament above 
kT = 0.5 keV, assuming an isobaric cooling flow model, is 
M = (2.46 + 0.13) x 10~ 2 M Q yr 1 . This cooling rate could 
only be reconciled with the apparent lack of kT < 0.5 keV 
plasma if the cold gas in the filaments absorbs the soft X-rays, 
modifying the observed spectral shape in a way that results in 
an apparent temperature cut-off. The X-ray spectral lines pro- 
viding the most diagnostic power at these temperatures are the 
Fe XVII and the O VII lines. The ionic fraction of Fe XVII 
is above 0.3 across a relatively large temperature range of 
0.2-0.6 keV, but the O VII line becomes observable only at 
kT < 0.4 keV. Absorption is expected to strongly suppress 
the O VII line emitted by the coldest gas phases. The Fe XVII 
lines will be affected by absorption to a much lesser extent. 
Therefore when modeling a spectrum that is affected by in- 
trinsic absorption, we may conclude that the cooling stopped 
at the temperature where we expect to start to see the O VII 
lin^] In Sect. |3] we show that the required amount of intrinsic 
absorption is N H = (1.24 + 0.12) x 10 21 cnr 2 . If the hot ICM 
indeed cools radiatively onto the cold filaments (the cooling 
takes place between the threads), then the large area covering 
fraction of the many small strands may provide a sufficient 
absorption column to the cooling X-ray emitting plasma to 
explain the apparent temperature floor at ~0.5 keV. 

Soft X-ray photons from the cooling X-ray plasma ab- 
sorbed by the cold gas would produce energetic photoelec- 
trons which would contribute to the heating and ionization of 
the cold gas (in a process similar to the secondary electron 
production by cosmic rays or impinging hot ICM particles, 
|Ferland et al.|20 09 ). Photoionization of emi ssion line nebu- 
lae by X-rays has be en considered before (e.g. Heckman et al. 
T989"l|Donahue & Voit|199T]|Crawford & Fabian|l9^2||Jaffe 
& Bremer|1997| ) and has been found insufficient to ionize and 
power the observed line emission. However, when consid- 
ering photoionization by primarily soft X-rays, from cooling 
plasma which surrounds the cold filamentary gas with a small 
volume filling fraction, the efficiency of this process increases 
significantly ( Oonk|20"TT i. According to our spectral model- 
ing, where we assumed that all the soft X-ray emission is due 
to radiative isobaric cooling, intrinsic absorption reduces the 
observed flux by ~ 3.6 x 10~ 14 erg s _1 cirT 2 . If the soft X- 
ray emitting plasma fills the space between the cold filaments 
and the cold gas absorbs approximately the same amount of 
energy in all directions, then the total X-ray energy absorbed 
by the cold gas is ~ 1.2 x 10 39 erg s _1 . With a total emit- 
ted energy from UV to sub-mm about 20 times larger than 
the observed Ho- luminosity of 3.0 x 10 38 erg s -1 , X-ray pho- 



: We note that only a small fraction of the observed X-ray emission orig- 
inates in the filaments and is affected by intrinsic absorption. The dominant 
emission component from the ambient ICM remains unaffected by cold gas 
in the filaments, and therefore fitting the spectra with a Galactic absorption 
as a free parameter will not result in a significantly higher best fit A?h- 
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toionization may, in principle, contribute a significant fraction 
of the energy needed to power the filaments. 

The cooling rates above were deduced under the assump- 
tion that all the soft X-ray line emitting gas originates from 
radiative cooling. The kT ~ 0.5 keV gas could, in principle, 
also be produced by mixing of the hot ICM with the cold fil- 
aments. More developed mixing might be responsible for the 
higher ratio of 0.5 keV thermal to Ha-+[N II] line emission 
in the southern filament (which is most likely falling back to- 
wards M 87), with respect to the northern filament (which 
is likely younger, currently being uplifted). Such mixing is, 
however, not expected to result in a gas with kT ~ 0.5 keV - it 
would result in a mixing l ayer with a wide range of tempera- 
tures ( Esquivel et al. 2006 ). Therefore, both in the cooling and 
in the mixing scenario for the origin of the soft X-ray emis- 
sion from the filaments, the explanation of the 0.5 keV cutoff 
requires intrinsic absorption by cold gas. 

The internal reddening in the optical band calculate d from 
the observed l ine ratio of Ha/HB = 4 in the filament ( |Ford &| 
Butcher|1979) is, E(B -V) = 0.28. Using the Galactic trans- 
formation between reddening and column density of neutral 
hydrogen, N H /E(B - V) = 5.8 x 10 21 cnr 2 mag" 1 (Bohlin 
et al. 1978[ ), we infer an intrinsic neutral hydrogen column 
density of Nu — 1.6 x 10 21 cuT 2 - sufficient to absorb the soft 
X-ray emission from gas with kT < 0.5 keV. We note that 
substantial reddening has been observed in a number of He? 
nebulae in brightest clust er galaxies of c ooling core clusters 
(e.g. |Allen et al.|1995||Crawford et al.p999] l. 

A similar temperature floor at ~0.5 keV is observed in the 
Perseus Cluster, Centaurus Cluster, 2A 03 35+096, Abell 262, 
Abell 3581, HCG 62, and Abell 2052 (|Sanders & F abian 
[2007i |Sanders et al.l[2008l [20091 |2010l |de Plaa et al.||2010| i 
and in the cooling uplifted multi-phase gas in the core of the 
galaxy cluster Sersic 159-03 (Werner et al. 201 1 1. The coolest 
~0.5 keV gas phases are co- spatial with filamentary optical 
emission-line nebulae in all of these systems. In the Perseus 
and Centaurus clusters, copious amounts of FIR line emis- 
sion has also been discovered, indicatin g the presence o f large 
reservoirs of cold atomic gas (Mitt al et al.|201 1 20 1 2) > . Fol- 
lowing the arguments outlined above, the low temperature X- 
ray cutoffs observed in those systems may therefore also, at 
least partly, result from intrinsic absorption. 

5. CONCLUSIONS 

We performed a multi-wavelength study of the emission- 
line nebulae located southeast of the nucleus of M 87. Our 
main conclusions may be summarized as follows: 

• We detect FIR [C II] line emission at 158 //m with 
Herschel PACS. The line emission is extended and co- 
spatial with optical Ha+[N II], FUV C IV lines, and 
soft X-ray emission. The observed [C II] flux implies 
1-5 x 10 6 M Q of cold gas within the field of view of Her- 
schel, of which ~ 0.4-2 x 10 6 M Q reside in the bright 
southern filament. 

• The filamentary nebulae contain multi-phase material 
spanning a temperature range of at least 5 orders of 



magnitude, from ~ 100 K to ~ 10 7 K. This material 
has most likely been uplifted by the AGN from the cen- 
ter of M 87. 

The thermal pressure of the 10 4 K phase appears to be 
significantly lower than that of the surrounding hot ICM 
indicating the presence of additional turbulent and mag- 
netic pressure in the filaments. If the turbulence in the 
filaments is subsonic then the magnetic field strength 
required to balance the pressure of the surrounding ICM 
is B ~ 30 - 70 fxG. 

The spectral properties of the soft X-ray emission from 
the filaments indicate that it is due to thermal plasma 
with kT ~ 0.5-1 keV, which is cooling by mixing with 
the cold gas and/or radiatively. Charge exchange can be 
ruled out as a significant source of soft X-rays. 

Both cooling and mixing scenarios predict gas with a 
range of temperatures. This is at first glance inconsis- 
tent with the apparent lack of X-ray emitting gas with 
kT < 0.5. However, we show that the missing very soft 
X-ray emission could be absorbed by the cold gas in the 
filaments with an average absorption column density of 
/Vh ~ 10 21 cm~ 2 , providing a natural explanation for 
the apparent temperature floor to the X-ray emission at 
kT ~ 0.5 keV. The internal reddening observed in the 
optical band indicates a similar level of intrinsic absorp- 
tion. 

The FIR through UV line emission is most likely pri- 
marily powered by the ICM particles penetrating the 
cold gas following a shearing induced mixing process. 
An additional source of energy may, in principle, be 
provided by X-ray photoionization from cooling X-ray 
emitting plasma. 

The relatively small line ratio of [O I]/[C II] < 7.2 indi- 
cates a large optical depth in the FIR lines. The large 
optical depth at FIR and the intrinsic absorption in- 
ferred from the X-ray and optical data all imply signif- 
icant reservoirs of cold atomic and molecular gas dis- 
tributed in filaments with small volume filling fraction, 
but large area covering factor. 
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